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■ We present results from lattice calculations on the thermodynamics of QCD at non- 

zero temperature and baryon chemical potential and discuss the role of resonances for the 
£\J \ occurrence of the transition to the quark-gluon plasma in hot and dense matter. Properties 

of a hadronic resonance gas are compared to lattice results on the equation of state at zero 
as well as non-zero baryon chemical potential. Furthermore, it is shown that the quark mass 
dependence of the transition temperature can be understood in terms of lines of constant 
energy density in a resonance gas. 
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§1. Introduction 

Today properties of strongly interacting matter at high temperature and non- 
zero baryon number density are analyzed within the field theoretic framework given 
by the theory of strong interactions - Quantum Chromo Dynamics (QCD). Two 
aspects of this theory, which describes the interaction among their elementary con- 
stituents, quarks and gluons, are of central importance for our understanding of 
the different phases of QCD at non-zero temperature and density as well as for our 
understanding of the experimentally observed spectrum of hadronic bound states - 
confinement and chiral symmetry breaking. The SUn(rif) x SU^nf) chiral symme- 
try, which is spontaneously broken in QCD with nj massless quark flavors, as well 
as the Z(N C ) center symmetry, realized in the gauge field sector of iV c -color QCD in 
the limit of infinitely heavy quarks, put strong constraints on qualitative aspects of 
the QCD phase diagram. 1 )' 2 ) In fact, quite general considerations concerning uni- 
versal properties of the QCD phase transition rely on these symmetries. Many of 
the emerging predictions concerning the order of the QCD transition and concerning 
universal properties in the vicinity of second order phase transitions in QCD have 
been verified over the last 20 years in numerical studies of lattice regularized QCD. 3 ) 
In our attempt to reach a deeper understanding of the physics behind the oc- 
currence of the QCD phase transition, mechanisms like the dual Higgs mechanism, 
monopole condensation or vortex percolation have been identified, 4 ) which char- 
acterize the drastic modifications of the QCD vacuum that occur at the critical 
temperature. However, none of these considerations provides insight into the ques- 
tion which properties of QCD set the scale and control quantitative aspects of the 
transition from hadronic matter at low temperature and density to the quark gluon 
plasma (QGP) at high temperature and/or density. 

Long before lattice calculations provided first evidence for a phase transition in 
strongly interacting matter 5 ) the inevitable need for critical behavior in hadronic 
matter has been discussed in the framework of resonance gas models. 6 ) It has been 
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Fig. 1. Sketch of the QCD phase diagram in the T-[ib plane (left) and the freeze-out curve (right) 
determined from a comparison of experimentally observed particle abundances to the abun- 
dances of hadrons in a hadronic resonance gas at temperature T and a baryon chemical potential 
Ms- 

noticed that ordinary hadronic matter cannot persist at arbitrary high temperatures 
and densities; 6 ) the copious production of resonances will lead to a natural end of the 
temperature and density regime in which hadrons can exist. It also has been sug- 
gested that the properties of dense hadronic matter and the role of resonances could 
be studied experimentally in heavy ion collisions. 7 ) Indeed, in such experiments it 
has been found that the abundances of various particle species are well described 
by a hadronic resonance gas model 8 ) which is characterized by two equilibrium pa- 
rameters, temperature, T, and baryon chemical potential, ^b- Both depend on the 
center of mass energy in these collisions; with increasing energy the temperature of 
the resonance gas increases and the relevant baryon chemical potential decreases. At 
RHIC the baryon chemical potential is quite small (/xb — 29 MeV) while the tem- 
perature of the resonance gas reaches 177 MeV. 9 ) This is in good agreement with the 
transition temperature to the quark-gluon plasma phase, T c = (173 ± 8 ± sys) MeV, 
found in lattice calculations 10 ) of 2-fiavor QCD at vanishing baryon chemical po- 
tential fiB- Recent exploratory lattice calculations also show that the transition 
temperature drops slowly with increasing baryon chemical potential 11 ) -13 ) and that 
the smooth transition at small \xb turns into a second order transition at a critical 
point 11 ) (T c , fj, B ) where T c is about 10% smaller than the transition temperature at 
[ib = and estimates for [i c B range from [i c B ~ 725 MeV 11 ) to [i c B ~ 420 MeV. 14 ) 
The anticipated QCD phase diagram and the freeze-out curve determined from the 
particle abundances observed in various heavy ion experiments are shown in Fiq. 1. 

In this paper we want to discuss evidence provided by lattice calculations for 
the contribution of hadron resonances to the thermodynamics of QCD. We will 
analyze the equation of state and the quark mass dependence of the QCD transition 
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Fig. 2. The transition temperature (left) as a function of the lightest pseudo-scalar meson mass 
(mps) in 2 and 3 flavor QCD and the energy density (right) in the SU(3) gauge theory (n/ = 0), 
as well as 2 and 3 flavor QCD for quark masses corresponding to mps - 770 MeV. 

temperature. Furthermore, we will address the question to what extend the freeze 
out temperature observed in heavy ion experiments is related to the phase boundary 
for the transition to the QGP calculated in lattice QCD. 

§2. Lattice results on T c and the equation of state 

2.1. (j, B = 

Studies of the quark mass dependence of the transition to the high temperature 
phase of QCD show that the transition temperature decreases gradually with decreas- 
ing quark mass*). Over a wide range of pion mass values, 300 MeV < mps < 2 GeV, 
the transition temperature depends linearly on mps and the slope seems to be more 
or less independent of the number of flavors (Fig. 2(left)). In Ref. 10 the mass 
dependence of T c has been parametrized as, 

T c (m PS ) = T c (0) + 0.04(l)m P5 . (2-1) 

For smaller pion masses one may expect that chiral symmetry leads to modifications 
of this linear relation* *); for mps > 2 GeV the transition temperature will approach 
a constant as all meson and baryon masses become larger than the (almost) quark 
mass independent glueball-masses and thus will decouple from the thermodynamics. 

The change in transition temperature between the light quark mass and infinite 
quark mass regime goes along with a change in e c /T^, the critical energy density 

*' Quark masses used in lattice calculations are bare parameters of the QCD Lagrangian which 
need to be renormalized. In order to avoid any discussion of the renormalization of quark masses it 
is more appropriate to discuss the quark mass dependence of thermodynamic observables in terms of 
a well controlled physical parameter. E.g. we will use here the lightest pseudo-scalar meson (pion) 
mass, mps- To express this in physical units (MeV) we use zero temperature lattice calculations of 
the string tension and set the scale by using y/a — 420 MeV. 

**' In the case of a second order transition in 2-flavor QCD one expects to find for small values 
of the pseudo-scalar mass, T c ~ mpj S , where 2/(36 ~ 1.1 is a combination of critical exponents 
of the 3-d, 0(4)-model. 




Fig. 3. The expansion coefficients C2 and C4 calculated in 2-flavor QCD for several values of the 
temperature. As in the [is — case the simulations have been performed with quark masses 
corresponding to mps — 770 MeV on a lattice of size 16 3 x 4. 



expressed in units of the transition temperature, by more than an order of magnitude 
(Fig. 2(right)), i.e. e c /T c 4 ~ (6 ± 2) for m PS ~ 770 MeV and e c /T c 4 ~ (0.5 - 1) 
for mps = 00. This large change reflects the large difference in the number of 
degrees of freedom which control the high temperature ideal gas limit. However, it 
does not at all suggest that the critical energy density itself changes significantly. 
In fact, when taking into account the shift in T c for both cases it seems that the 
critical energy density itself does not change much. Although statistical errors are 
still large both cases are consistent with a critical energy density in the range e c ~ 
(0.5- 1.0) GeV/fm 3 . 

2.2. p B > 

Lattice simulations at non-zero baryon chemical potential*) generally suffer from 
the problem that the fermion contribution to the QCD partition function is no 
longer represented by a positive definite quantity; the fermion determinant becomes 
complex. This excludes the application of standard Monte Carlo simulation tech- 
niques. 15 ) This problem has, however, been avoided in recent studies of the de- 
pendence of the transition temperature on the chemical potential 11 ^ 13 - 1 as well as 
calculations of the equation of state 16 )' 17 ' at non-zero baryon chemical potential by 
using extrapolation techniques applied to numerical results obtained at \xb = 0. 

At fixed temperature and small values of the chemical potential the pressure 
may be expanded in a Taylor series around \i q = 0, 

1 00 

Yi = Vr^ z = ^ CniT) ' (2 ' 2) 

?1=0 



*' We use here the notion of baryon chemical potential [ib although lattice calculations are 
performed in terms of the quark chemical potential fi q = /j.b/3. 
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Fig. 4. The left hand part of the figure shows the change in pressure due to a non-zero quark 
chemical potential, fi q = /is/3 calculated in 2-flavor QCD using a Taylor expansion to order 
(Hq/T) 4 . As in the fi B = case the simulations have been performed with quark masses 
corresponding to mps — 770 MeV. The right hand part gives results on the /^s-dependence of 
the transition line together with the chemical freeze-out line of Fig. 1 which corresponds to a 
line of constant energy per particle in a a hadronic resonance gas. 



where the expansion coefficients are given in terms of derivatives of the logarithm 

1 d n In Z 

of the QCD partition function, c n (T) = — ttw — , - . The series is an even series 

n. o{fiq/ 1 ) 

in (fiq/T). The first coefficient, cq, just gives the pressure studied for some time in 
finite temperature lattice calculations at fi q = 0. The first non-zero coefficient, c%, 
is proportional to the quark number susceptibility at [i q = which also has been 
studied in the past. 18 ) The relevant expansion coefficients entering the calculation of 
the pressure and quark number susceptibility in a Taylor expansion of the partition 
function up to 0(fj lq ), 17 ' 1 

Ap(T,Hq) _ p{T, fi q ) - p(T, 0) _ ( Hq\ 2 , //V 4 



%L = 2c 2 + 12cM) 2 , (2-3) 



rp2 Z \ T 

are shown in Fig. 3. In Fig. 4(left) we show corresponding results on the change 
in pressure, Similar results have been obtained using a reweighting technique. A 
comparison with the pressure calculated at [i q = shows that at n q /T = 1 and 
for T > T c the enforced presence of a non-zero baryon number adds about 30% to 
the overall pressure in the system. This also confirms that the Taylor expansion 
converges rapidly up to fi q /T ~ 1. 

The reweighting approach 11 )' 12 ) as well as the analytic continuation of results 
obtained in simulations with an imaginary chemical potential 13 ^ have also been used 
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to determine the /^-dependence of the transition temperature. These results are 
shown in Fig. 4(right). We note that for small values of the chemical potential 
the lines shown in this figure do not correspond to a phase transition but rather 
characterize a rapid but smooth crossover from the hadronic to the plasma phase. 
Also shown in the figure is an estimate of the chiral critical point, n ) i.e. a second 
order phase transition point, (T c , n° B ) ~ (160 Me V, 725 MeV), in the QCD phase 
diagram. For fis > /J.% the transition is expected to become a first order phase 
transition. 

Although the different methods used to determine T c (/zb) do seem to give com- 
patible results, a detailed quantitative comparison is difficult as all calculations per- 
formed so far have used different discretization schemes and/or quark mass values. 
Moreover, the estimates have partly been obtained on rather small lattices with 
unimproved gauge and fermion actions and/or too large quark masses. The current 
quantitative results thus need to be improved and confirmed in future calculations. A 
first analysis of the quark mass dependence of the transition line for 3-flavor QCD 19 ) 
seems to indicate that the //^-dependence becomes stronger with decreasing quark 
mass and that the chiral critical point in the physically realized case of (2+l)-flavor 
QCD shifts to smaller values of the baryon chemical potential, 14 ) \i c B ~ 420 MeV. 

§3. Thermodynamics of the hadronic resonance gas 

3.1. fi B = 

In Hagedorn's approach to the thermodynamics of strongly interacting matter 6 ) 
critical behavior arises because any increase in energy of an ensemble of strongly 
interacting hadrons is predominantly used to generate new resonances rather than 
transforming it into kinetic energy of the constituents. Energy thus is not used for 
heating up the system any further. The exponentially rising spectrum of resonances, 
p(m) ~ exp (bm), leads to the occurrence of a critical (limiting) temperature, T c res — 

1/6- 

A similar mechanism is known to lead to critical behavior in purely gluonic 
systems, i.e. in the SU(N C ) gauge theories. Here the fluctuations of color flux 
tubes (string) lead to an exponentially rising excitation spectrum which again leads 
to critical behavior. Calculations of the critical temperature, T CiStr i ng , within the 
Nambu-Goto model yield, 20 ) 

^string _ / 3 . . 

which only depends on the space-time dimension, d, and, in particular, is indepen- 
dent of the color degrees of freedom. Lattice calculations of the phase transition 
temperature of SU(N C ) gauge theories in 3 and 4 space-time dimensions, indeed, 
yield critical temperatures, which are in good agreement with the string model pre- 
dictions. Some results from lattice calculations 21 ) are summarized in Table I and 
compared with the string model prediction. This clearly suggests that resonances 
play an essential role in triggering the occurrence of the deconfinement transition, 
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Table I. Lattice results on the deconfinement temperature in units of the square root of the string 
tension of SU(N C ) gauge theories in d space-time dimensions. The last row gives results obtained 
from a large-d analysis of string models. 20 ' 
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although the order of the transition and, in those cases where the transition is second 
order, also universal properties of the transition are controlled by the global Z(N C ) 
center symmetry of the SU(N C ) gauge theories. It thus is interesting to explore what 
role the (exponentially) rising hadronic resonance spectrum plays for the occurrence 
of the transition to the plasma phase in the physical, light quark mass regime. 

A first hint at the importance of resonances for the occurrence of the QCD 
transition may be obtained from the energy density at T c . In the limit of nf massless 
quark flavors the sector of massless Goldstone bosons would contribute to the energy 
density with e/T c 4 = (n 2 — l)7r 2 /30, if this contribution can be approximated by an 
ideal gas of non-interacting bosons. The calculations of the energy density shown in 
Fig. 2 (right) have been performed with quark masses which correspond to a pion mass 
of about 770 MeV. Their contribution to the energy density thus is exponentially 
suppressed. A free gas of massive relativistic particles would contribute 
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3T K 2 {km/T) 
km 



+ Ki{km/T) 



(3-2) 



to the energy density. Here n = —1 for bosons, +1 for fermions and g is the degen- 
eracy factor of the particle state. For m/T c ~ 4 this yields e/T^ ~ 0.083(n 2 — 1) 
which shows that quite a few hadronic degrees of freedom are needed to saturate the 
value of the energy density found in lattice calculations at T c . 

In a gas of non-interacting resonances energy density and pressure are given as 
sum over the single particle contributions, e.g. e = X^ el ( m «)> where i labels the 
masses, nii, of experimentally known mesons and baryons. In order to compare the 
resonance gas model with lattice results we, however, have to take into account that 
the latter are not obtained from calculations performed with the physically realized 
light quark mass spectrum. A corresponding analysis has been performed in Refs. 
22,23. The resulting comparison between a modified resonance gas model and lattice 
data for (2+l)-flavor QCD 10 ) is shown in Fig. 5. Similar agreement has been obtained 
for 2-flavor QCD by suppressing the contribution of "strange" hadrons. This shows 
that resonance can account for the rapid rise of the energy density observed in lattice 
calculations for T<T C . 
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Fig. 5. The energy density in (2+l)-flavor QCD 
(e - 3p)/T 4 , (right) as a function of T/T c . 
resonance gas with hadron masses adjusted to 
been performed with too heavy up and down 
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(left) and the deviation from ideal gas behavior, 
The solid line gives the result obtained from a 
compare with the lattice calculations which have 
quarks. 22 ' 



3.2. hb > 

The resonance gas model does make quite stringent predictions for the ther- 
modynamics at non- vanishing chemical potential. 22 - 1 ' 23 ' 1 As all baryons in QCD are 
heavy compared with the temperature regime of interest, i.e. mbaryon/^ ^ 5 for 
T < T c , the contribution of baryons to the thermodynamics can be handled in the 
Boltzmann approximation. This leads to a factorization of the temperature and 
fugacity, exp(fj,s/T), dependence of thermodynamic observables. The pressure of a 
gas of baryons and their resonances thus can be written as 

PbF^b) =Fb{t) cosh(/Wr) ; (3 . 3) 

where Fb{T) is given by a sum over all baryons and their resonances, 

9i /m.-\ 2 



i 

As the mesonic part of a gas of non-interacting resonances does not depend on the 
baryon chemical potential the change in pressure due to a non-vanishing chemical 
potential is entirely determined by the baryonic sector, Ap = p(T, /^b) — p(T,0) = 
Pb{T,[ib) — Pb(T,0). This, of course, also holds for derivatives with respect to \xb 
at fixed temperature. The factorization of the T and (//e/T^-dependent terms then 
leads to simple relations among various thermodynamic observables. For instance, 
one finds for the baryonic susceptibility, 

XB = ( d 2 p(T,fi B ) \ 

T 2 U/WT) 2 J Tflxed 

= ^ (1 - cosh" WT))- 1 • (3-5) 

This relation is particularly interesting as it suggests that the ratio Ap/T 2 x q = 
Ap/9T 2 XB does not depend on details of the hadron mass spectrum and thus can 
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Fig. 6. The ratio of pressure and quark number susceptibility \ q = 9xb versus temperature for 
fixed values of the quark chemical potential fJ, q /T — {ib/3T (left). The horizontal lines are 
the results of hadron resonance gas model calculations. 23 ' The points are the lattice data from 
Ref. 17. While the dashed-dotted curves represent the complete expression given in Eq. 3-5 the 
dashed curves give the result of a Taylor expansion performed to the same order as that used in 
the lattice calculations. The right hand part of the figure shows the quark number susceptibility 
versus T/T c for different values of the quark chemical potential. The lines give results from the 
resonance model calculation expanded in a Taylor series and truncated at the same order as 
used in the lattice calculation. 



directly be compared with lattice calculations performed with vm-physically large 
quark masses. Such a comparison is shown in Fig. 6(left). The agreement between 
the lattice results and the resonance gas model relation is quite striking, although it 
has to be noted that the current lattice analysis has been performed using a Taylor 
series expansion to order only. This truncation, of course, eliminates any true 
density driven singular behaviour and, for instance, cannot lead to a divergent sus- 
ceptibility. Large fluctuations which could occur in the vicinity of a phase transition 
are suppressed. This also is true for the resonance gas model approach, which does 
not lead to true singular behaviour as long as it is formulated in terms of a finite set 
of resonances. 

The resonance gas model also yields a reasonably good description of the tem- 
perature dependence of thermodynamic quantities at fixed fugacity. This, of course, 
involves information on the explicit form of the resonance spectrum. Taking again 
into account that lattice simulations so far have been performed with unphysical 
quark mass values one can directly compare the quark number susceptibilities calcu- 
lated in the resonance gas model and on the lattice. This is shown in Fig. 6(right). 
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Fig. 7. Transition temperatures in 2 and 3 flavor QCD versus the lightest pseudo-scalar meson 
mass. Both quantities are expressed in units of the square root of the string tension. The solid 
lines correspond to temperatures at which the energy density of a 2 and 3 flavor resonance gas 
attains the value e = 0.8/GeV/fm 3 . The former consist of hadrons without strangeness content 
while the latter also includes strange mesons and baryons. 

§4. The critical temperature 

The good agreement between lattice results on the QCD equation of state in the 
low temperature phase and the resonance gas model raises the question whether we 
can also understand the dependence of the QCD transition temperature on the quark 
mass in terms of a resonance gas. The similarity of the "critical" energy densities 
found in the infinite quark mass limit and for (moderately) light quark masses sug- 
gests that the transition to the quark gluon plasma phase occurs at approximately 
constant energy density, irrespective of the quark mass values or the physical masses 
of the hadronic resonances. This assumption is the basis for the comparison of tran- 
sition temperatures in 2 and 3 flavor QCD with lines of constant energy density 
calculated in a resonance gas model. 22 ) Results are shown in Fig. 7. Up to pseudo- 
scalar masses mps — (3 — 4)y / <7 ~ (1.3 — 1.7) GeV the agreement with the hadronic 
resonance gas is quite reasonable. For larger pseudo-scalar masses the glueball sec- 
tor does start to play an increasingly important role as the heavy hadrons decouple 
from the thermodynamics and the "lighter" glueballs yield the largest contribution 
to the energy density. Aside from the inclusion of glueball states in the resonance 
gas 22 ) one eventually also has to take into account that thermal effects may strongly 
influence the glueball spectrum close to the transition temperature. 24 ) Doing so the 
qualitative features of the quark mass dependence of T c can be modeled in the entire 
mass regime. 

Finally, let us consider the lines of constant energy density in the T-[i B plane. 
Separating the meson and baryon contribution to the energy density, 

e(T,n B ) = e M (T) + e B {T)cosh(n B /T) , (4-1) 

where c_b(T) is the baryonic contribution to the energy density at fj, B = 0, one now 
can follow the strategy applied also in the lattice calculations and determine the 
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lines of constant energy density from a leading order Taylor expansion. Expanding 
around the transition point at vanishing chemical potential, (T c , [xb = 0), yields, 



We note that the denominator on the right hand side is just the specific heat in a 
resonance gas at (T c , hb = 0). This shows that one does expect a quite weak depen- 
dence of the transition temperature on the chemical potential, if the transition at 
[iB = is close to a second order transition. Although this is in qualitative agreement 
with the weak /xe-dependence of the transition line observed in lattice calculations 
it seems that the resonance gas still leads to a somewhat stronger variation with \x b ■ 
A more direct comparison in the limit of physical quark masses would be desirable 
here. 



We have shown that the hadronic resonance gas model is able to describe quite 
a few quantitative and qualitative results obtained in lattice calculations on the 
thermodynamics of the low temperature hadronic phase of QCD as well as basic 
properties of the transition line to the high temperature phase. This suggests that 
the copious production of hadronic resonance indeed plays an important role in 
triggering the transition to the quark gluon plasma phase of QCD. Unfortunately, 
a direct comparison of lattice calculations with an non-interacting gas of resonances 
given directly in terms of the experimentally known hadron spectrum is at present 
not yet possible. If one wants to compare the current lattice calculations, which still 
are being performed with too heavy quark masses, with the hadronic resonance gas 
some phenomenological input on the quark mass dependence of hadron resonances 
is required. However, it will soon be possible to further reduce or even eliminate 
the present ambiguities. With the availability of a new generation of Teraflops 
computers for lattice QCD it soon will become possible to perform studies of QCD 
thermodynamics with an "almost" realistic spectrum of quark masses. 



This paper summarizes talks given at the conferences on "Finite Density QCD" 
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would like to thank the organizers, in particular A. Nakamura and H. Suganuma, for 
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